Long-range effects on superdiffusive solitons in anharmonic chains 
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Studies on thermal difFusion of lattice solitons in Fermi-Pasta-Ulam (FPU)-like lattices were re- 
cently generalized to the case of dispersive long-range interactions (LRI) of the Kac-Baker form. 
The position variance of the soliton shows a stronger than linear time-dependence (superdiffusion) 
as found earlier for lattice solitons on FPU chains with nearest neighbour interactions (NNI). In con- 
trast to the NNI case where the position variance at moderate soliton velocities has a considerable 
linear time-dependence (normal diffusion), the solitons with LRI are dominated by a superdiffusive 
mechanism where the position variance mainly depends quadratic and cubic on time. Since the 
superdiffusion seems to be generic for nontopological solitons, we want to illuminate the role of the 
soliton shape on the superdiffusive mechanism. Therefore, we concentrate on a FPU-like lattice 
with a certain class of power-law long-range interactions where the solitons have algebraic tails in- 
stead of exponential tails in the case of FPU-type interactions (with or without Kac-Baker LRI) . A 
collective variable (CV) approach in the continuum approximation of the system leads to stochastic 
integro-differential equations which can be reduced to Langevin-type equations for the CV position 
and width. We are able to derive an analytical result for the soliton diffusion which agrees well with 
the simulations of the discrete system. Despite of structurally similar Langevin systems for the two 
soliton types, the algebraic solitons reach the superdiffusive long-time limit with a characteristic t^'^ 
time-dependence much faster than exponential solitons. The soliton shape determines the diffusion 
constant in the long-time limit that is approximately a factor of tt smaller for algebraic solitons. 



INTRODUCTION 



Nonlinear excitations like solitons and discrete 
breathers have been drawing great attention over recent 
years. They present very robust solutions of nonlinear 
partial differential equations and nonlinear lattice mod- 
els which are often used to describe a rather broad set of 
physical systems 0, 0, ■ Nonlinear lattices like Klein- 
Gordon (KG), Fermi-Pasta-Ulam or the discrete nonlin- 
ear Schroedinger (DNLS) lattice (and their counterparts 
in the continuum approximation) present all to often a 
strong approximation of the system of interest and it is 
sometimes not clear to what extend solitons or discrete 
breathers are also relevant in more realistic models. In 
the case of biomolecules, the aspiration to highlight their 
basic functionalities with computer simulations is a chal- 
lenging task for the computational science nowadays. If 
we aim to understand principally the role of nonlinear ex- 
citations in biomolecules, we have to investigate simpler 
models than the numerical ab initio calculations. Often 
used in this context is the Davydov model for energy 
and charge transport in proteins or the Peyrard-Bishop- 
Dauxois (PBD) model which describes the melting and 
the denaturation of DNA 0, 0- There are convincing 
evidences that some features of nonlinear excitations in 
this simplified models are relevant in explaining the func- 
tionality of biomolecules. Recently, pump-probe mea- 
surements showed that the lifetime of NH stretching 
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bands in the model protein ACN is (with about 18 ps) in 
good agreement with numerical calculations of the Davy- 
dov model The observation that the PBD model 
succeeds in determining the thermally induced openings 
of the DNA strand at functionally relevant sites for the 
DNA transcription is also an example for the rele- 
vance of nonlinear lattice models in biology. 
Nowadays, many physicists who work in this field try 
to extend and improve the basic models in order to 
make them more realistic. The spatial structure of 
biomolecules, thermal fluctuations, damping and long- 
range effects stemming from Coulomb or dipole-dipole 
interactions have certainly a great influence on the non- 
linear excitations and the statistics of the system. It 
is known that long-range effects can change the fea- 
tures of solitons and discrete breathers qualitatively when 
the LRI exceed some critical value [H [l^. In FPU- 
like chains and nonlinear Schrodinger models, they give 
rise to new types of solitons which can coexist at the 
same value of the spectral parameter |0. Controlled 
switching between such soliton states was recognized as 
a possible mechanism for energy transport and storage in 
biomolecules In the past, nonlinear KG lattices with 
LRI were frequently studied in investigations of a num- 
ber of physical phenomena such as dislocations in solids, 
charge density waves, absorbed layers of atoms or domain 
walls in ferromagnets and ferroelectrics (see references in 
Elllllil)- The effect of thermal fluctuations and LRI 
are usually regarded separately 0, 0, |23| because of 
their complexity. Nevertheless it is known that LRI can 
have very interesting effects on the thermodynamics of 
many different physical systems |^ . 
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The Kac-Baker form presents a spatially exponential cou- 
pling between different particles which is quite often cho- 
sen for the inclusion of LRI in nonlinear lattice models 
p^, '2I, '2^ . In the recent publication lattice solitons 
on FPU- like chain with Kac-Baker LRI in the presence of 
a thermal reservoir were investigated. In the continuum 
approximation of the system and with a collective vari- 
able (CV) approach, it was possible to derive a Langevin 
set for the soliton position, similar to the solitons on FPU 
chains [2^. The solitons on FPU chains (with or with- 
out Kac-Baker LRI) show superdiffusive behaviour. The 
position variance of the soliton shows besides the linear 
term in the time-dependence also quadratic and cubic 
terms. In it was demonstrated, that the diffusion 
mechanism for NNI and for additional LRI are quite sim- 
ilar because the same CV ansatz was used. In the case 
of Kac-Baker LRI, this approximation is valid unless the 
soliton velocities do not reach the critical velocity, where 
the soliton begins to develop a cusp 12]. 
In order to check the influence of the soliton shape on the 
diffusion of the soliton, we choose a power-law coupling. 
This leads in the continuum approximation (CA) to a 
Benjamin-Ono (BO) equation, which is known to posses 
algebraic soliton solutions |22|. The first studies in that 
direction date back to Ishimori [2^ who studied anhar- 
monic chains with Lennard- Jones (2n,n) intermolecular 
potential and showed that the dynamics is governed by 
the Benjamin-Ono equation in the case rt = 2 or by the 
Korteweg-de Vries equation for n > 4. In the case of 
cubic and quartic nearest neighbour interactions, a re- 
ductive perturbation method yields a KdV-|-mKdV and 
Benjamin-Ono equation wich was shown to possess exact 
nonsingular rational solutions |26|. 



The model covers the physical situation of dipole-dipole 
(s — 5) and Coulomb interactions (s — 3) between the 
particles on the chain, if we restrict ourselves to small 
relative displacements. The equation of motions in rela- 
tive displacement coordinates can be obtained from the 
Lagrangian of the system L = T — Unn — Ulr 



bS) = V'{w,+^)-2V'{w,) + V'{w,_^) 



J 



'Ji+m ) 



m/0 



m 



(5) 



When we want to consider the effect of a thermal bath, 
we can add damping and noise terms to the original equa- 
tion such that they fulfill the ffuctuation-dissipation 
theorem {D'^y = 2vhyT with ks set to unity) 



m{t) = V'{w,+i)~2V'iw,) + V'iw,^i) 

jsr^ im- Wi+m) , . . o • , • ^ 

-J} j r- h i^/iy(Wi+l - 2Wj + Wi_i) 



(6) 



We choose hydrodynamical damping which depends on 
the relative displacement velocities of the particles and 
presents an intrinsic damping mechanism of the system. 
The widely used Stokes damping is not appropriate for 
pulse solitons because one obtains an imaginary disper- 
sion relation the long- wave region of the Fourier spectrum 
what causes deformations of the soliton (2^. We go to 
the continuum limit [wilt) — > w{x,t) , ^i(t) C,ix,t) , 
f{x+m, t) e™^^/(a;, t)] and find the partial differential 
equation (PDF) 



II. THE MODEL 

Our model is a one-dimensional chain of equally spaced 
particles of mass m (to = 1) with an interatomic spacing 
a (a = 1). We denote the displacement of the particle 
n from its equilibrium position (absolute displace- 

ment coordinates) and the difference Wn — Wn+i — Wn as 
relative displacement coordinates. The potential consists 
of a part Unn similar to the original potential chosen by 
Fermi, Pasta and Ulam and a long-range potential Ulr 
with a power-law dependence of the harmonic coupling 
constant Jnm 



N 



Unn = ^ V{un+i{t) - Un{t)) 



Ulr^ 



2 \ m — n [ 



(1) 

(2) 

(3) 
(4) 



dtw{x, t) 



2{cosh{d^) - l)V'{w{x,t)) 
1 — cosh{mdx) 



2jj: 



w{x, t) + yhyd^dtw{x, t) 



•m—l 

+ ^WydU{x,t) . 



(7) 



The further treatment of the system depends on the value 
of s. In [23 it is shown that in the case s > 5 the CA 
yields a Boussinesq equation similar to the result with- 
out LRI but with a different dispersion parameter. For 
3 < s < 5 it was proved that the soliton tails are no 
longer exponential but algebraic and that for s < 3.5 
an energy gap between the soliton states and the plane 
wave spectra appears. For s — A, the equation of motion 
becomes a Hilbert-Boussinesq equation which can be re- 
duced to the integrable Benjamin-Ono form, which has 
algebraic soliton solutions. 

Since the soliton equation in |23| was of Boussinesq-type 
and covered the limit of nearest neighbour interactions, 
we do not expect fundamental new effects for s > 5. In 
order to investigate the behaviour of solitons with alge- 
braic tails in the presence of thermal ffuctuations, we 
restrict ourselves to the case s = 4 where the Hilbert- 
Boussinesq equation and the resulting Benjamin-Ono- 
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type solitons present a promising system to achieve ana- a ^ _ ^2 _ ^-ly/aS^, , J_a4 /-.qn 

lytical results. The perturbed Hilbcrt-Boussinesq equa- (^1^,0^;- ^ g ^l^xJ + 24 ^ ^ 



tion 



d^w ~ c^dlw - Xdtw - —1-L{dlw) + 25^^^ = 

= vhydldtu + VD^d.^ix, t) (8) 

with = (1 + j7rV6) und A = (1/12- J/24) foUows from 
expanding in the operators for the nearest-neighbour 
part and for the long-range part Q{s,dx) = 2^,^(1 — 
coshimdx)) / m'^ 

2(co.M5x)-l) = 2(^1 + 1 + ...) (9) 



where H{f{x)) = l/irP f{y)/{y — x)dy denotes the 
Hilbert transform and P the Cauchy principal value. 
Similarly like in [2^, we will have to rewrite the soli- 
ton equation ^ in absolute displacement coordinates in 
order to find its Lagrangian density (in the case Vhy = 0). 
Notice that the A-term was neglected in reference [T5| . 



FIG. 1: (Color online) Soliton solutions for power-law (s = 4) 
and Kac-Baker long-range interactions. The two solitons in 
the left and the right panel each have the same velocity Co 
and the same energy H. 




III. COLLECTIVE VARIABLES 

Jtt 

The Hilbert-Boussinesq equation in absolute displace- d^u - c^dlu - Xd^u —H{dlu) + 2dxudxxU 

ment coordinates u[x,t) 
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= Uhyd^dtu + VD^d^^ix, t) (11) 
can be derived for v^^ = from the Lagrangian density 



(12) 



As a CV ansatz with X{t) and cr(t), we will use the 
Benjamin-Ono-type soliton shape 



u{x, t) = ~ AoArctan[a{t){x - X{t))] 
W(X, t) = — - 



(13) 



i + a{t)^x-x{t)y ■ 

One can check that this ansatz yields the correct values 
Ao = ^+ AXao , a, = -^{v^ - c') (14) 

3 TTJ 

(compare with [ISj) when one minimizes the action for 



L 



I 



dxC 



c A Tra TT 



4 

^ + y/D^y I dx 



- -A^a^ + 
4 8 

Act 

1 + a^x- Xit))' 



A^m^a'^ 

48 



(with A = and V D^v = 0) for the soliton parameters 
(7, A and X(i) assuming a coherent excitation with 
constant velocity {& = A = X = v). 

To include the damping, we proceed like in with 
the generalized Hamilton principle of Ostrovsky et al. 



S<C> d<C> d d<C> 



5X, 



dx, at dx. 



= (16) 



where $ = v^^u^xt is a dissipative field and the brackets 
signify a spatial integration over x. 

In Fig. n we show two solitons with Co — 1.01 (left panel) 
and Co = 1.03 (right panel) on FPU chains with power- 
law LRI (s = 4, J = 0.7877, c = 1.51516) and for the 
Kac-Baker LRI (described in l24j). The values of a and 
J for the solitons with the exponential tails in the Kac- 
Baker case were chosen to produce a solution with the 
same velocity and soliton energy H . The power-law LRI 
yield solutions with lorentzian shape (algebraic solitons) 
whereas the Kac-Baker LRI yield solitons which are well 
approximated by a sech (exponential solitons). 



IV. LANGEVIN SYSTEM 

The eqs. (|16f) together with the ansatz (|13|l leads to 
following stochastic equations for the collective variables 



a and X : 



aX = 



4VD'"^f7^ 



i>2 



with 



Ao^ 
AoO Acr^ 



dx- 



x^{x,t) 



(1 + Cr2s2)2 



ttA„ 



dx 



^hy V, 3 
——A cr 



(1-^(72^2)2 



(17) 



Ar, 



Jtt 



-t- 4Acro , X 



X{t) 



(18) 



In the previous step, we neglected small terms ^ , ^ X 
and 7 in order to obtain technically less comprehensive 
stochastic equations. The justification of this step was 
checked by numerical simulations of the complete system 
and was already successfully applied in the case of Kac- 
Baker LRI H. 

We want to rewrite the second equation as an equation 
for X. We expand the appearing square root for the small 
value V D^y to first order and use that the velocity of the 
damped soliton 



(19) 



which can be approximated by w, the start velocity of the 
soliton (because c < Vd < v, c ^ v). 

We end up in a system of stochastic integro-differential 
equations for the CV: 



Al 
A2 



dx 



Bn 
B 



22 



with 



Al 



Bu 



A2 = Vd 

{x-X{t)) 



B. 



ttAov (l + a2(a;-X(t))2)2 
(x-X(i))2 



4/0^0-2 



ttAov {l + a^{x-X{t)Yy 



(20) 

(21) 
(22) 

(23) 



We proceed like in |33| in order to find a statisti- 
cally equivalent Langevin system to H20|l (with the same 
Fokker-Planck equation in the Stratonovich interpreta- 
tion |33l |) which is more convenient for further numerical 
and analytical studies. The Langevin-system with two 
independent Gaussian white noise processes reads: 



with 



ai 



bii 
622 



(24) 



AA^nv- 



0.2 = Vd 



hi = — ; — 1=— , 022 



Ar, 



/TTV 



y^Aoy/nv 



(25) 
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The corresponding system in |2J| depends similarly on 
the soliton parameters inverse width, Ac, velocity v and 
on the temperature VZ3- Only numerical constants and 
the velocity Vd are different. The broadening of the soli- 
ton due to the damping follows the same manner as for 
exponential solitons in j2^l25l|. The result H24|) manifests 
the statement made in |2j|, that the LRI and the result- 
ing soliton shape determine the velocity Vd{<j{t)) of the 
damped soliton. 

V. SMALL-NOISE EXPANSION 

We proceed with a small-noise expansion for the pa- 
rameter V D^y like in E^. l25l | in order to derive an ana- 



lytical expression for the soliton diffusion. The Langevin 
system (124(1 is formally very similar to the sytem in [23 | 
for Kac-Baker LRI, only numerical constants and the soli- 
ton parameters are different. For algebraic solitons, we 
can approximate H19|l with Vd — c + Aq/^c a{t) because 
of the smallness of A. The analytical solution in this case 
is technically less extensive and leads practically to the 
same results. The approximation for Vd is not possible 
in the case of Kac-Baker LRI in because the long- 
range forces yield a larger term which is quadratic in the 

inverse soliton width. 

The result in the first order of VL)^ reads (see Ap- 
pendix) 



Var[X(t)] = D 



Ag[-16(-1 + V^) + (7? - 1)(18 - 10^?; + (7? - 1)(7? + 2))] 2c2z.^^7?(~l + V^-(7?-l)(7y + l)) ^| 

(26) 



with 



VI. SIMULATIONS 



T] = l + tr = l + (rlvhyt . (27) 

The time-dependence of the position variance depends 
on the time scale tr, which explains that the large time 
limit (and the superdiffusive behaviour) sets in earlier for 
high- velocity solitons (large values of (Tq) for the same 
damping constant. The corresponding result for Kac- 
Baker LRI is even more complex, which stems from the 
longer expression for Vd 24]. The time scale which de- 
scribes the broadening of the exponential soliton and the 
diffusion is tr = LGi^hyJo where j{t) is the inverse width 
of the exponential soliton 

We are interested in the soliton diffusion for large times 
tr- Substituting H27|) into (|26|) and looking for the leading 
order of tr yields: 

Jjhy 



Var[Xit)],„f 



-t 



1.5 



(28) 



hy 



which differs from the result in |24j only by the numerical 
constant. 

For very small times, the result (|26|l gives a linear time- 
dependence 

J^hy 



Var[X{t)]o = -- 



rt 



(29) 



The results of the small-noise expansion shows the same 
characteristic dependencies on the soliton parameters as 
for Kac-Baker LRI in [2^ • For intermediate times, the re- 
sult 1(261) describes the transition when the position vari- 
ance turns from the small-time dependence into the long- 
time limit wherefore a stronger increase in time (mainly 
~ t^) can be observed. 



The simulation of system ll^ was performed in the 
same manner as described in |24| but with power-law 
long-range coupling instead of the Kac-Baker LRI. The 
position variance of the soliton was calculated from 100 
different realizations of the chain. We used the soli- 
tons depicted in Fig. as initial conditions of the 
chain. The soliton with Co = 1.01 (ctq ~ 0.1) is much 
broader than the lattice spacing whereas the soliton with 
Co = 1.03 (cto ~ 0.3) is rather discrete in the soliton cen- 
ter. The coupling parameter is always fixed at the value 
J = 0.7877 which yields the sound velocity c = 1.51516, 
the same value was used in [2^ . In order to compare the 
results of the Kac-Baker LRI with the power-law LRI 
and to be able to draw meaningful conclusions , we also 
use the same damping constant Vhy = 0.01 and the same 
temperature T = 0.0001 as in Ref. 0. 

The analytical result for the position variance l(26() 
agrees rather well with the simulation results (Fig. O. In 
general, one expects the analytical results to be slightly 
smaller than the simulation results because the possible 
contribution of phonons is not regarded in the calcula- 
tions of Section III. It was shown for a related short-range 
model that the phonons have indeed a perceptible influ- 
ence on the soliton diffusion "sF]. The analytical result 
for the low- velocity soliton diffusion for Co — 1.01 is too 
small for small times. Similar discrepancies appeared for 
Kac-Baker LRI in and they were interpreted as the 
influence of phonons which cause higher values in the 
simulations. The analytical result for the high-velocity 
soliton Co = 1.03 deviates from the simulation results for 
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FIG. 2: (Color online) Position variance of a low-velocity 
(1.01) and a high-velocity (1.03) soliton on a anharmonic 
chain with power-law long-range interactions (J = 0.7877, 
s = 4). We compare the simulation results with equation 
12till . the result of a small-noise expansion of the Langevin 
system 11241 of the CV theory. 



1=0.787679, s=4, co=1.01, T=0.0001, v=0.01 




1=0.787679, s=4, co=1.03, T=0.0001, v=0.01 



Simulation 
YaT[X(t)] "small noise expansion" 



responding time for Cq = 1.01 (uo « 0.1) is t > 10000. 
Therefore, we will check the validity of the long-time ap- 
proximation result of the position variance for the soliton 
with Co — 1.03. In deriving the long-time limit 128|) of the 
position variance, we recognized that the approximation 



Var[X{t)]u 



D 



ti 



(30) 



describes the simulation results very well for all times t. 
Neglecting the y^-term for large tr yields immediately 
the result p^ . 



FIG. 3: (Color online) Position variance of a high-velocity 
soliton (co — 1.03) on an anharmonic chain with power- law 
long-range interactions (J = 0.7877, s — 4). The simu- 
lation and two different long-time approximations for I26L 
Var[X{t)],r,f ^ and Var[X{t)]h JSOJ, are shown. For 
t — 10000 {tr ~ 9), the long-time limit is not yet reached, but 
the simulation results approach this function for long times. 



1=0.787679, s=4, c =1.03, T=0.0001, v=0.01 




4000 6000 8000 10000 



The comparison between the simulation results and 
Var[X{t)]fi or Var[X{t)]inf in Fig. O shows that there 
is practically no linear dependence on t similarly as for 
high-energy solitons with exponential tails with Kac- 
Baker LRI in 24]. The superdiffusion with stronger- 
than-linear terms in the time-dependence of the position 
variance is the dominating mechanism for high-energy 
sion. When the solitons are very narrow, the width of the solitons. The time-dependence oi Var[X{t)]h describes 
soliton changes quite fast and the first-order correction the simulation results very well for all times and the soli- 
of a{t) can fail. Such deviations were also observed for ton can be described by a time-dependence according to 
FPU chains with or without Kac-Baker LRI (24„,3£|. (|28|l (~ f'^/^) for large times. For lower soliton ener- 

The analytical result for the soliton diffusion depends gies like for Co = 1.01 (Fig. ^ the normal contribu- 
only on the time scale = vty'^'tt and for large times tions on yar[X(t)] are stronger and a comparison with 
the diffusion should approach the result 1)28(1 from be- Var[X(t)]inf yields a worse agreement because in this 
low. Since the time scale tr depends on the width of case the limit of large tr can not be used for t < 10000. 
the soliton, the long-time approximation for the position In the limit of short times {tr << 1), the quality of the 
variance for high- velocity solitons will be valid for smaller result for the linear contribution in time H29|l is very simi- 
values of t than for low-velocity solitons. For Cq = 1.03 lar to the Kac-Baker case and yields values for I^ar[X(i)] 
(cto ~ 0.3) tr > 1 is valid for t > 1111 whereas the cor- which are smaller than the simulation results. For the 




large times. This feature seems to stem from our restric- 
tion to first order perturbations in the small-noise expan- 
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soliton with Co = 1.01 the slope of (|29() is approximately 
half of the value one would get from a linear fit of the 
simulation results for small times. 



VII. COMPARISON 

Up to now, we have mentioned the similarities between 
the thermal diffusion of solitons with algebraic and ex- 
ponential tails on FPU-like chains. Now, we want to 
demonstrate the main difference, namely the different 
time-dependence of the superdiffusion in the two cases. 
We want to directly compare the results for yar[X(t)] 
and their time-dependence for two solitons with the same 
energy and velocity (co = 1.03) but one on a chain with 
power-law LRI and the other on a chain with Kac-Baker 
LRI. For the two solitons, normal diffusion contributions 
are negligible and the best fits in Fig. ^ demonstrate 
that the position variance of algebraic solitons increases 
with ~ i"^/^ whereas for the sec/i-shaped solitons, the in- 
crease goes with a quadratic (~ t'^) and a cubic (~ t^) 
term. The small-noise expansion results for both soliton 
types have the same long-time behaviour (~ t^^"^) (com- 
pare with |23|) which is reached by the algebraic soliton 
in Fig. 21 for t > 8000, whereas the exponential soliton 
is for times t < 10000 still in the regime where the long- 
time limit is not applicable. 

The Langevin systems for both soliton types are identi- 
cal apart from the numerical constants and the values of 
the soliton parameters like the width. For Kac-Baker 
LRI, the inverse soliton width 70 = 0.132 at i = 
is distinctly smaller than for the power-law case with 
(To — 0.279 for the same energy and soliton velocity. Since 
the time scale tr depends quadratically on the inverse 
width, it is not surprising that the long-time limit for ex- 
ponential solitons is valid after much longer times. This 
argument is decisive for different solitons of the same 
type. It follows directly from the small-noise expansion 
that the long-time limit for the algebraic soliton with 
Co = 1.01 appears for a^{co = 1.03)/al{co = 1.01) w 5.44 
longer times than for the algebraic soliton with Co = 1.03. 
In order to estimate when the exponential soliton reaches 
the long-time asymptotics, we have to concentrate our 
studies on the discussion of the result yar[X(i)] from the 
small-noise expansion for algebraic and exponential soli- 
tons (Fig. |SJ|. Simulations for times much longer than 10** 
would be too time-consuming. The double-logarithmic 
scales in Fig. O make clear that the Kac-Baker soli- 
ton reaches the long-time limit Var[X(t)]inf after much 
longer times than the power-law soliton. The time win- 
dow where a stronger time-dependence than ~ t^^^ ap- 
pears is approximately ten times larger for the Kac-Baker 
soliton. 

The influence of the soliton shape determines the diffu- 
sion constant in Var[X{t)]inf , which is approximately a 
factor of TT larger for exponential solitons. One can equiv- 
alently state that the longer transient times to the long- 
time limit for the exponential soliton is a consequence of 



FIG. 4: (Color online) Simulation results for the position vari- 
ance of two solitons, one on an anharmonic chain with power- 
law long-range interactions (J = 0.7877, s = 4) and the other 
soliton on an anharmonic chain with Kac-Baker long-range 
interactions for T = 0.001. The interaction radius and the 
coupling for the Kac-Baker case were chosen in such a way 
that the two solitons have the same velocity Co = 1.03 and 
the same energy. In order to demonstrate the difference in the 
time-dependence for the two solitons, we added a best fit for 
the two simulation results for Var[X{t)]. One can clearly see 
that the position variance of the algebraic soliton approaches 
a ~ t^'^ time-dependence, whereas in the Kac-Baker case, 
the time-dependence can be fitted by a quadratic and a cubic 
term. 



H=1.09, c^=1.03, c=1.51516, T=0.0001, v=0.01 




2000 4000 6000 8000 10000 



the higher value of the diffusion constant. This remark- 
able difference for the two soliton types results only from 
the different soliton shapes because the numerical con- 
stant in H28|l is determined by the numerical constants in 
the Langevin system H24|l which follow from spatial inte- 
grations over the soliton profile in the CV procedure. 
The time evolution of the velocity of the energetically 
equivalent algebraic and exponential soliton, Vd{(T{t)) and 
Vdilit)), is rather different despite of the same start 
value and the structural similarities for Vd- In Fig. (BJ 
we present the simulation results for the mean velocity 
{X{t)) of the two solitons in comparison with the zeroth 
order expressions for Vd from the small-noise expansion: 

Vd{a^°Ht)) e.nd Vdh^^'Ht)) M- 

These hold for the zero temperature case but describe 
the mean velocity of the two noisy solitons already quite 
well. We can directly see that the power-law soliton ap- 
proaches the long-time limit with small soliton velocities 
earlier than the Kac-Baker soliton. The slower dynamics 
of the Kac-Baker soliton is caused by its broad, exponen- 
tial shape and the resulting slower time scale tr which 
controls the time-evolution of "f^°\t). The dependence 
of Vd on a term ~ 7^ for exponential solitons (which 
is negligible for algebraic solitons) influences the soliton 
velocity (and therefore the superdiffusion), especially for 
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FIG. 5: (Color online) Var[X{t)] and Var[X{t)]^r,f from a 
small-noise expansion for the two solitons of Fig. 0] The 
algebraic soliton reaches the long-time limit Var[X{t)]i„f 
much sooner than the exponential soliton with the same ini- 
tial energy H = 1.09 and velocity Co = 1.03. The result 
Var[X{t)]inf shows the same time-dependence but the diffu- 
sion constant for algebraic solitons is approximately a factor 
of TT smaller. 



FIG. 6; (Color online) The simulation results for the mean 
soliton velocity (X) in case of a noisy chain with T — 0.0001 
for energetically equivalent algebraic and exponential solitons 
with Co = 1.03. The zeroth-order results (without noise) for 
the soliton velocities Vd{u^'^\t)) and Vd{'y^'^\t)) agree rather 
well with the simulation results and demonstrate the different 
dynamics of the two soliton types. 
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H=l.ro,c =1.03, c= 1.51516 



100 




VarlXMl, 




^, Kac-Baker: 1=0.505541, a=0.765 






VarlXMl, 




^, power-law: 3=4, 1=0.787679 






^"IXMlj 


BCac-Bater 


1=0.505541, a=0.765 






_- V>rIX(t)l_^ 


power-law 


3=4, 1=0.787979 




<dX/dt> Kac-Baker: J=0.505541, 0=0.765 
o <dX/dt>, power-law: s=4, J=0. 787679 

v/Y'°'(t)) 
--- vW°'(t)) 




pronounced LRI character and for small times when the 
long-time limit is not yet reached. 



VIII. CONCLUSIONS 

We investigated the diffusion of algebraic solitons on 
FPU-like chains with power-law long-range interactions 
and derived analytically a result for the position vari- 
ance which agrees well with the simulations of the dis- 
crete system. The results clearly reveal that the soli- 
ton diffusion and the damping-induced reshaping of the 
soliton depends on the time scale tr = a'^Vhyt which is 
given by the soliton width and the damping constant. 
For very broad solitons, the damping effects are quite 
weak because the particles in the soliton profile have 
moderate relative velocities and evolves slowly in t. A 
small-noise expansion of the Langevin system yields an 
analytical result for the soliton diffusion which depends 
on tr and predicts a superdiffusive long-time behaviour 
~ t^/^. This result makes clear why the superdiffusion, 
which sets in for large tr, dominates the diffusion of high- 
energy solitons (cr~^ small). These basic features were 
also found in earlier studies for solitons with exponential 
shape. The superdiffusion of lattice solitons can therefore 
be regarded as generic when the lifetime of the solitons is 
much larger than l/a^Vhy- On chains with NNI, this sit- 
uation can often not be reached for low-energy solitons 
because they are rapidly destroyed by the fluctuations 
and the damping. The LRI stabilize the soliton and lead 



to higher soliton energies and lifetimes which makes the 
superdiffusion the dominating mechanism. 
The shape of the soliton does not change the basics of 
the diffusion process. The Langevin systems for the col- 
lective variabels (inverse width and soliton position) are 
technically similar, only numerical constants and the de- 
pendence of the soliton velocity on the inverse width 
are different. We proved in simulations and by analyz- 
ing the analytical result that despite of the similarities, 
the superdiffusive behaviour of the two soliton types ex- 
hibits striking differences even when the soliton energy 
and velocity coincide. The algebraic solitons are gen- 
erally narrower than exponential solitons which leads 
to faster dynamics for the time scale tr and an ear- 
lier validity of the long-time limit with the characteristic 
~ t^/^ time-dependence. The broadness of the exponen- 
tial solitons, the corresponding slower dynamics for tr 
and the stronger dependence on higher orders of the in- 
verse width in the soliton velocity shift the ~ t'^/^ region 
to larger times t which leads to the observed character- 
istic quadratic and cubic time-dependence in the simu- 
lations. The perhaps most remarkable effect of the dif- 
ferent soliton shapes is that the exponential solitons show 
larger position variances than their algebraic counterparts 
for sufficiently long times. The diffusion constant in the 
long-time limit is for algebraic solitons approximately a 
factor of TT smaller than for exponential solitons. This 
effect is only caused by the different soliton shapes which 
determine numerical constants in terms which describe 
the coupling between the soliton velocity and the fluctu- 
ations of the soliton width. Further studies for different 
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values of s would be helpful to explain the influence of 
power-law interactions on the soliton diffusion. Espe- 
cially the diffusion of solitons on chains with power-law 
LRI for the case s < 3.5 could be interesting because a 
gap between the plane-wave spectra and the soliton en- 
ergies appears. 



APPENDIX: SMALL-NOISE EXPANSION 

The small-noise expansion of the Langevin set l|24|l fol- 
lows reference Q and was already successfully applied 
to the problem of the diffusion of low-velocity solitons 
in the FPU system without long-range interactions [25|. 
We seek an asymptotic solution of the form 



ea 



(1) 



(A.l) 
(A.2) 



where e is a small parameter which is formally introduced 
to consider the influence of the noise terms as small per- 
turbations {V D^y ^ e). The contribution of the drift 
term in the different orders of e are calculated following 
the rule 



(A.6) 



These equations describe the damped soliton. The damp- 
ing induced broadening of the soliton is formally equiva- 
lent to the case of exponential soliton solutions in ,241 • 
the expression for the soliton velocity it is obvious that 
the damped soliton gets broader (smaller a) which leads 
to a slowdown of the soliton. In the limit of long times, 
the soliton width diverges and the soliton approaches the 
velocity of sound. 

The first order corrections due to the noise read: 



2 Ao^nv 



(A.7) 



4c V(T(o)^/^Ao^; 



If we substitute the result for ct^"-' = ctq/ \/v^^^a1t + 1 
into (|A.8p we calculate the first-order expression for the 
dislocation of the soliton due to the noise 



a,(a(t)) = a,(a(") + ^e'"a(")(t)) 

m— 1 

In order to minimize the technical efforts, it is advisable 
to approximate the velocity 02 = Vd by the following 
expression 



(A.4) 



Peril) 

Vd 



* 3 /•*' 3 

dtr]' 5 / dt"ri"^dWi{t") 
Jo 

rA i\'~'dW2{t') . (A.9) 

(A.IO) 



1 + UhyO-lt' 



After some straightforward calculations we can calculate 
the position varinace 



Var[Xi{t)] = \im)Xi{t)Xi{s){ (A.ll) 



which was checked to be appropriate for typical soliton 
parameters J, v. The equations in the order e'^ read: 



(A.5) 



which yields the result H26|l where we can finally set e — 
1 which is justified if the noise in the system is sufficiently 
small. 
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